Field locked to Fock state by quantum feedback with single photon corrections 
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Fock states with photon numbers n up to 7 are prepared on demand in a microwave super- 
conducting cavity by a quantum feedback procedure which reverses decoherence-induced quantum 
jumps. Circular Rydberg atoms are used as quantum non-demolition sensors or as single photon 
emitter/absorber actuators. The quantum nature of these actuators matches the correction of single- 
photon quantum jumps due to relaxation. The flexibility of this method is suited to the generation 
of arbitrary sequences of Fock states. 

PACS numbers: 42.50.Pq, 42.50.Dv, 03.67.Pp 



The preparation of non-classical field states and their 
protection against decoherence is an important aspect of 
quantum physics and of its application to information 
science. Among the proposed methods, including error 
correction [ij, decoherence- free subspaces m and reser- 
voir engineering [3|, quantum feedback [^-[g is particu- 
larly promising. Its principle is to drive the quantum 
system towards a target state by the repeated action of 
a sensor-controller-actuator loop. The sensor performs 
quantum measurements and provides information to the 
controller. Taking into account the back-action of the 
measurement, the controller then estimates the system's 
state and programs the actuator to drive the system as 
close as possible to the target. A given feedback al- 
gorithm can in principle protect a wide class of target 
states. The operating point can be changed at any time 
and the system driven through a programmed trajectory 
in its Hilbert space. 

Photon number (Fock) states are appealing targets for 
quantum feedback operation. They combine a theoret- 
ical and intuitive simplicity (|n) is an eigenstate of the 
field Hamiltonian with n quanta) with intrinsically non- 
classical features (their Wigner functions take negative 
values for n > 1). When coupled to an environment, Fock 
states lose rapidly their non-classicality with a time con- 
stant Tn = Tc/n, where Tc is the lifetime of the field mean 
energy [7|. Hence, large n Fock states, whose decay time 
scale is much shorter than Tc, are important tools for 
the exploration of decoherence at the quantum/classical 
boundary (8|. 

The most intuitive algorithm for quantum feedback 
generation and protection of Fock states combines sin- 
gle photon emitters/absorbers together with a Quantum 
Non-Demolition (QND) photon counting sensor. Start- 
ing from vacuum, repeated photon emissions make the 
field climb the ladder of Fock states until the sensor rec- 
ognizes that the target is reached. When an environment- 
induced single-photon quantum jump occurs, the sensor 
detects it and the controller triggers the action of either 
an emitter or absorber to correct for it. 



FIG. 1. Scheme of the experimental set-up. 



Resonant actuators and QND sensors have already 
been used separately to prepare Fock states in a time 
short compared to T„. Repeated photon emission by 
resonant actuators has been achieved in cavity- [9| and 
circuit-QED experiments [10|. QND photon counting 
by repeated atomic sensor measurements has collapsed 
a coherent field into a Fock state with random n and 
revealed its subsequent quantum jumps 
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Single- 
photon-actuators and QND sensors have never been com- 
bined in a feedback loop. The only steady-state quantum 
feedback procedure so far has stabilized Fock states with 
a coherent source as actuator [12]. The performance of 
this scheme was limited by the inability of the classical 
actuator to correct for a quantum jump in a single step. 
We report here the implementation of a quantum feed- 
back algorithm for the preparation and protection of Fock 
states up to n = 7 in a microwave cavity QED experi- 
ment using circular Rydber g at oms as QND sensors and 
as single-photon actuators [13|. The experimental set- 
up is depicted in Fig. [T] The microwave field (frequency 
v — b\ GHz) is stored in the superconducting Fabry- 
Perot cavity C {Tc = 65 ms) cooled at 0.8 K (average 
number of blackbody photons nth = 0.05) [ij]. Rydberg 
atoms are prepared in B into the circular state \g) of prin- 
cipal quantum number 50 by a pulsed process repeated 



at Ta = 82 fxs time intervals. This laser- and r.f.-induced 
process also selects the atomic velocity {v = 250 m/s). 
The small number of atoms in each sample obeys a Pois- 
son law. 

A voltage V applied across the cavity mirrors tunes 
the atoms in or out of resonance with C by Stark-shifting 
the transition between \g) and the nearest higher- lying 
circular state |e). After crossing C, the detection of each 
atom in |e) or \g) by state-selective field- ionization in D 
provides binary information. 

The cavity C is placed between two low-Q cavities Ri 
and R2 used to manipulate the atomic state by classical 
resonant microwave pulses (sources Si and S2). The ex- 
periment is controlled by a computer K (ADwin Pro-II), 
which processes binary information from D and uses it 
to adjust for each sample the settings of Si, S2 and V. 
When K decides to send a sensor in C, V is set to make 
the atom- field interaction dispersive and Si and S2 are set 
for 7r/2 pulses. The R1-R2 combination is then a Ramsey 
interferometer providing QND information on n by mea- 
suring the dispersive phase shift experienced by the atom 
in C 111. When K decides to send an emitting actuator, 
it applies with Si a 7r-pulse realizing the \g) — >■ |e) trans- 
formation in Ri and it sets V so that the atom interacts 
resonantly with C during an adjustable time. The source 
S2 is switched off and the atom is directly detected in 
D. For an absorbing actuator, Si and S2 are switched off 
and the resonant interaction time is controlled by V. 

The first task of K is to estimate the field state af- 
ter detection of each sample, based on all available in- 
formation. Since the initial vacuum and the actuators 
bear no phase information, the field density matrix re- 
mains diagonal in the {|7i)} basis and K needs only 
to update the photon number distribution p{n) after 
each detection. From p(n), K evaluates the distance 
d — J2n{n — ntYp{n) = An^ -f- (n — ntY to the tar- 
get Fock state \nt) (An^ and n are the photon number 
variance and mean value, respectively). 

After detection of a sensor in state j = or 1 (|e) or 
\g), respectively), p{n) becomes, according to Bayes' law 
and within a normalization factor, p{n) x ■Ks[j\n), where 
■ns{j\n) is the conditional probability for detecting the 
sensor in j when there are n photons. Ideally, 7rs(j|n) = 
[l-|-cos($o?^ + ¥'r— J7'')]/2, where $0 is the phase shift per 
photon accumulated by the atomic coherence in C and ipr 
the adjustable phase of the Ramsey interferometer |12| . 
In the experiment, due to various imperfections, 7rs(j|n) 
becomes [1 + jbs + Cg cos(<I>o'^ + Vi- ~ J''')]/2, where the 
offset bs = 0.02 ±0.002 and the contrast Cs = 0.75 ±0.03 
are determined in calibration experiments. 

The resonant actuators and C perform a Rabi os- 
cillation between the joint states \e,n) and \g,n + 1) 
at a frequency proportional to \/n + 1. After detec- 
tion of an actuator, prepared in j and detected in k 
{j,k = 0,1 for |e) and \g), respectively), the pho- 
ton number distribution p{n) becomes, within a nor- 



malization, p{n + j ~ k) X 7ra(j, fcjri + j — k), where 
Tra{j,k\n + j — k) is the conditional probability of the 
j ^ k transition in the field oin + j — k photons. Ideally, 
■^aU^Mn+j-k) = {l+cos[notj^/n - fc + l+(j-/c)7r]}/2. 
Here, rio/27r ~ 47.9 kHz is the vacuum Rabi oscilla- 
tion for an atom at cavity center; tj — tg or tg where 
te and tg are effective interaction times taking into ac- 
count the atomic motion through the Gaussian structure 
of the mode in C. Experimental imperfections affect the 
contrast of these oscillations in an n-dependent way and 
slightly offset them. The actual 7ra(j, A;|n) are obtained 
by recording Rabi oscillations in calibration experiments. 

The state estimation must also take into account the 
measured detection efficiency {rid = 0.25). For instance, 
if no atom is detected, it is possible that the sample was 
actually empty [p{n) being then unchanged] or that it 
contained one, or even two undetected atoms. The con- 
troller K updates p{n) as a weighted average of the mod- 
ified photon number distributions corresponding to all 
possible events leading to the actual detection. The prob- 
abilities of these events are inferred by Bayes' law, know- 
ing rjd and the mean atom number m in a sample [15|. 
Undetected sensors have no effect on state estimation, 
since a QND atom does not modify p{n) on the average. 
Undetected actuators must be considered, the modified 
p{n) being a weighted average of the ones which would 
be obtained if the missed atom would have been detected 
in |e) or \g). If two actuators are crossing C together (0, 
1 or 2 being detected), the updating formulas take into 
account the modified Rabi oscillation and the possibility 
of simultaneous two-photon emission or absorption. 

The state estimation includes cavity damping during 
the time interval Ta towards the thermal equilibrium at 
the mirrors temperature [15|]. Finally, K takes into ac- 
count the feedback loop delay, i. e. the effect of the yet 
undetected actuator samples that are on their way be- 
tween C and D, by tracing over their states. 

Getting information on p{n) requires many sensors 
(each providing one bit), whereas adding or subtracting a 
photon ideally demands a single atom. Thus, a feedback 
loop consists of a sequence of Ng sensor samples followed 
by Nc < Ns control samples whose operating mode is 
decided by K. 

After each detection and state estimation, K makes 
or updates a decision about the fate of all control sam- 
ples located immediately before Ri or between Ri and C. 
In the first case, K decides upon the sample mode (sen- 
sor, emitter or absorber). For actuators between Ri and 
C, K uses additional information acquired after selection 
of their modes (emitter/absorber) to decide to let them 
interact resonantly with C as previously planned or to 
cancel their action by keeping them off-resonance. The 
controller estimates the average photon number distribu- 
tion resulting from the pending interactions of the control 
samples with C for each of these choices and selects that 
corresponding to the minimum distance d to the target. 
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FIG. 2. Single realization of the feedback experiment with 
nt = 4. The frames present versus time, from top to bottom, 
the detected sensor states (upwards bars for e, downwards 
bars for g), the distance d to the target, the actuators sent 
by K (red bars for emitters, blue bars for absorbers) and the 
photon number distribution p(n) inferred by K (color/gray 
scale) together with its average value (solid black line). 



The free experimental parameters are the phase-shift 
per photon $0: the Ramsey phase (pr, the actuator inter- 
action times tj, the number of sensor and control sam- 
ples Ns and Nc and the corresponding average num- 
ber of atoms per sample rus and rric. The phase-shift 
$0 — 0.252 X TT rad, corresponding to an atom-cavity de- 
tuning 5/2tt = 244 kHz, is set close to 7r/4 allowing K to 
distinguish among eight different photon numbers llj. 
The Ramsey interferometer phase ipr — 7''/2 — $o"-t is set 
by fine Stark-tuning of the atomic frequency. It corre- 
sponds ideally to TTs{j\nt) = 1/2 and provides the best 



sensitivity to photon number measurements around nt- 
We chose irig ~ 1.3 and rric = 0.5 (the lower value 
for the control samples reduces the probability of two- 
photon emission/absorption). The other parameter val- 
ues, te = 1.67r/rJo\/"t + 1, tg = 2.47^/^io^/^, Ns = 12 
and Nc = 4 are optimized by numerical simulations. 

The value of tg is close to 27r/i7o\/"-t + 1- This cor- 
responds to a "trapping state" condition [16|, for which 
an emitting actuator would ideally leave the target state 
invariant. Due to the finite contrast of the experimental 
Rabi oscillations, the emission probability does not cancel 
at the trapping state condition. Choosing a slightly lower 
value for t^ maintains a relatively small unwanted emis- 
sion probability for n = nt while optimizing the proba- 
bility of correcting emissions when n = nt — 1. Similar 
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FIG. 3. Probabilities of the choices made by K for the con- 
trol atoms (emitter: dashed red line; sensor: solid green line; 
absorber: dash-dot blue line) as a function of the estimated 
mean photon number n (data inferred from 4000 realizations 
of the experiment over 140 ms with nt —4). 



arguments explain qualitatively the value of tg, slightly 
larger than that corresponding to a trapping state con- 
dition. 

Figure[2]shows the data of a single realization of the ex- 
periment with nj = 4. It presents, as a function of time, 
the detected sensor states, the estimated distance d, the 
controller decisions to send emitter or absorber actuator 
samples, and finally the evolution of the photon num- 
ber distribution estimated by K together with its average 
value. Starting from vacuum at i = 0, emitting samples 
are repeatedly sent until d comes close to zero. The pho- 
ton number distribution is then peaked on ti = n^, with 
p{nt) ~ 0.8 — 0.9. Around i = 50 ms, a downwards quan- 
tum jump to 71 = 3 triggers the sending of few emitter 
samples, which rapidly restore the target state. Close to 
t = 70 ms, another downwards jump is over-corrected, 
leading to n = 5. Absorbers are then sent until restora- 
tion of the target. Four thousand similar trajectories 
have been recorded for each value of nt from 1 to 7. 

In order to get an intuitive insight into the workings of 
the feedback, we plot for n^ = 4 in Fig. [3] the fractions of 
emitters, absorbers and sensors chosen by K in the control 
samples, versus the mean photon number n estimated at 
a given time. The mere inspection of this figure leads to 
a simple rule. When n < nt — 0.4 (resp. n > nt + 0.6), 
K essentially programs emitter (resp. absorber) samples 
and for nt — 0.4 < n < nt + 0.6 it rather decides to 
send sensor samples, which do not affect n on average, 
but contribute to reducing An and thus d. The domain 
of n values in which sensors are preferred has a 5n sa 1 
width and is centered on a value slightly larger than nt, 
reflecting the fact that field relaxation alone reduces n. 
The fractions of actuators and sensors vary rapidly at 
the boundary of these domains, in narrow ranges of n 
values. Our feedback algorithm thus essentially operates 
according to the simple intuitive procedure outlined in 
the introduction. When it estimates that n falls outside 
the Sn = 1 range around nt, it concludes that a quantum 
jump is likely to have happened and attempts to correct 
for it by adding or subtracting a photon. 

The performance of the feedback procedure is obtained 




FIG. 4. Histograms of photon number distribution as a function of the target photon number nt- (a) Reference Poisson 
distribution with rit photons on the average, (b) Photon number distribution p(n) measured by an independent QND process 
after interrupting the feedback loop at t = 140 ms. (c) Photon number distribution p(n) measured when K estimates that 
p{nt) > 0.8. For (b) and (c), p(n) is measured from no to no + 7, with no = for nt < 5, no = 2 for nt > 5. 



by reconstructing the average p{n) distribution, indepen- 
dently from the estimation made by K, and by comparing 
it to the reference Poisson distribution with nt photons 
on the average shown in Fig.|4l[a). The loop is interrupted 
at i = 140 ms. We then send a few QND sensor samples 
with Ramsey phase settings optimizing the discrimina- 
tion of 8 photon numbers ^llj . From 4000 realizations of 
this experiment, we reconstruct p(n), plotted in Fig.|4ljb) 
between no and nn_+ 7 around nt, by a maximum like- 
lihood procedure [8|. The measured p{nt), i. e. the fi- 
delities with respect to the target state, are about twice 
those of the corresponding Poisson distribution and the 
p(n) distributions are clearly sub-Poissonian. 

A better performance is obtained when using the state 
estimation by K for interrupting the feedback at a proper 
time. The histograms in Fig. Sl^c) present for each rit the 
p(n) distributions obtained with 4000 sequences inter- 
rupted when K estimates that p(nt) > 0.8. The obtained 
fidelity is close to the expected 0.8 value up to nt = 3 
and remains always larger than that of the histograms 
in Fig. \Mjo). From a statistical analysis of these trajec- 
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FIG. 5. Programmed sequence of Fock states. The target nt 
varies stepwise as indicated by the thin line. The p(n) dis- 
tribution inferred by K is shown in color/gray scale, together 
with its average value (thick line). 



tories, we also infer the convergence time towards the 
target. For nt — 3, 63% of the trajectories have reached 
the 0.8 fidelity threshold within 27 ms. This time is twice 
as short as that observed with classical actuators |12| . 

To illustrate the fiexibility of the procedure, K is pro- 
grammed to change nt in real time according to a preset 
sequence, {3, 1, 4, 2, 6, 2, 5} for the single trajectory pre- 
sented in Fig. [5] The controller switches from one nt 
value to the next when it estimates that p{nt) > 0.8. It 
accordingly adapts the Ramsey phase fr to the new tar- 
get. The p{n) distribution and its average n (thick black 
line) follow rapidly the changes in the target (thin blue 
line) . 

The operation of the feedback mechanism can also be 
viewed as that of a quantum micromaser [16[ pumped by 
the actuators and actively locked to a Fock state. Based 
on information provided by the sensors, the controller 
adjusts in real time the atomic medium gain by choos- 
ing the fraction of emitters and absorbers. In ref. [1^ the 
controller was rather injecting coherent microwave pulses 
into the cavity. The mismatch between the classicality of 
the source and the quantumness of the target limited the 
procedure to low photon number states. Here in contrast, 
the quantum nature of the source allows us to correct 
more rapidly and precisely decoherence-induced quantum 
jumps. We lock efficiently the micromaser operation to 
higher Fock states, making them available for fundamen- 
tal studies and quantum information experiments. 
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